Convective coupled map for simulating spatiotemporal chaos in flows

A coupled map lattices with convective nonlinearity or, for short, Convective Coupled Map (CCM) is proposed in this paper to simulate spatiotemporal chaos in fluid flows. It is found that the parameter region of spatiotemporal chaos can be determined by the maximal Liapunov exponent of its complexity time series. This simple model implies a similar physical mechanism for turbulence such that the route to spatiotemporal chaos in fluid flows can be envisaged. The study is supported by “Nonlinear Sciences Project” from the State Science and Technology Commission of China.     

Crack in a material-bond lattice

A square-cell lattice is considered consisting of point masses at its knots connected by linearly elastic bonds of nonzero density. Steady-state crack propagation is studied. A general relation between the knot mass and the bond mass is assumed; however, a detailed analytical examination is made for the material-bond lattice with no concentrated masses. It is assumed that the crack divides the bond in half, and the broken bonds remain in the lattice structure. In this model, the fracture energy of the bond is ignored, and hence the local fracture energy of the lattice is zero. The classical formulation in terms of critical stresses is accepted. The macrolevel energy release does exist. The macrolevel energy release rate as a function of the crack speed is found and compared with that for the massless-bond lattice of the same averaged density. While in the main, the dependencies for these two models are similar, there are some essential differences. For the lattice with no concentrated masses this function appears discontinuous. There exists a region where the crack speed is insensitive to the variation of the macrolevel energy release rate. The admissible regions of the crack speeds for the considered two lattice models differ greatly. For the massless-bond lattice this region is rather wide, while for the other it is very narrow. Mathematically, it is of interest that some details of the factorization depend on whether the ratio of the crack speed to the wave speed is rational and, if so, whether it can be represented as a ratio of two odd numbers.     

Numerical Investigations of the Steady State Relative Permeability of a Simplified Porous Medium

The purpose of this paper is to investigate, by flow simulations in a uniform pore-space geometry, how the co and countercurrent steady state relative permeabilities depend on the following parameters: phase saturation, wettability, driving force and viscosity ratio. The main results are as follows: (i) with few exceptions, relative permeabilities are convex functions of saturation; (ii) the cocurrent relative permeabilities increase while the countercurrent ones decrease with the driving force; (iii) with one exception, phase 2 relative permeabilities decrease and phase 1 relative permeabilities increase with the viscosity ratio M=₁/₂; (iv) the countercurrent relative permeabilities are always less than the cocurrent ones, the difference being partly due to the opposing effect of the viscous coupling, and partly to different levels of capillary forces; (v) the pore-level saturation distribution, and hence the size of the viscous coupling, can be very different between the cocurrent and the countercurrent cases so that it is in general incorrect to estimate the full mobility tensor from cocurrent and countercurrent steady state experiments, as suggested by Bentsen and Manai (1993).(Now at AS Norske Shell, Norway.) e-mail:     

Partial Uniqueness and Obstruction to Uniqueness in Inverse Problems for Anisotropic Elastic Media

We consider the inverse problem of identifying the density and elastic moduli for three-dimensional anisotropic elastic bodies, given displacement and traction measurements made at their surface. These surface measurements are modelled by the dynamic Dirichlet-to-Neumann map on a finite time interval. For linear or nonlinear anisotropic hyperelastic bodies we show that the displacement-to-traction surface measurements do not change when the density and elasticity tensor in the interior are transformed tensorially by a change of coordinates fixing the surface of the body to first order. Our main tool, a new approach in inverse problems for elastic media, is the representation of the equations of motion in a covariant form (following Marsden and Hughes, 1983) that preserves the underlying physics.In the case of classical linear elastodynamics we then investigate how the type of anisotropy changes under coordinate transformations. That is, we analyze the orbits of general linear, anisotropic elasticity tensors under the action by pull-back of diffeomorphisms that fix the surface of the elastic body to first order, and derive a pointwise characterization of parts of the orbits under this action. For example, we show that the orbit of isotropic elastic media, at any point in the body, consists of some transversely isotropic and some orthotropic elastic media. We then derive the first uniqueness result in the inverse problem for anisotropic media using surface displacement-traction data: uniqueness of three elastic moduli for tensors in the orbit of isotropic elasticity tensors. ^†Partially supported by an MSRI Postdoctoral Fellowship. Research at MSRI is supported in part by NSF grant DMS-9850361. This work was conducted while the first author was a Gibbs Instructor at Yale University. ^‡Partially supported by an MSRI Postdoctoral Fellowship, and by NSF grant DMS-9801664 (9996350).     

Extending Circle Mappings to the Annulus

We show that any monotone degree one C mapping of the circle can be realized as the induced mapping on the boundary by a C area preserving diffeomorphism of the open annulus. The circle mapping may have critical points, allowing the possibility of Denjoy behavior on the boundary.     

扁球面网壳的混沌运动研究

在圆形三向网架非线性动力学基本方程的基础上,用拟壳法给出了圆底扁球面三向网壳的非线性动力学基本方程．在固定边界条件下,引入了异于等厚度壳的无量纲量,对基本方程和边界条件进行无量纲化,通过Galerkin作用得到了一个含二次、三次的非线性动力学方程．为求Melnikov函数,对一类非线性动力系统的自由振动方程进行了求解,得到了此类问题的准确解．在无激励情况下,讨论了稳定性问题．在外激励情况下,通过求Melnikov函数,给出了可能发生混沌运动的条件．通过数字仿真绘出了平面相图,证实了混沌运动的存在．     

Studies on Bifurcation and Chaos of a String-Beam Coupled System with Two Degrees-of-Freedom

In this paper, research on nonlinear dynamic behavior of a string-beam coupled system subjected to parametric and external excitations is presented. The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system. The Galerkin's method is employed to simplify the governing equations to a set of ordinary differential equations with two degrees-of-freedom. The case of 1:2 internal resonance between the modes of the beam and string, principal parametric resonance for the beam, and primary resonance for the string is considered. The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system. Based on the averaged equation obtained here, the techniques of phase portrait, waveform, and Poincare map are applied to analyze the periodic and chaotic motions. It is found from numerical simulations that there are obvious jumping phenomena in the resonant response–frequency curves. It is indicated from the phase portrait and Poincare map that period-4, period-2, and periodic solutions and chaotic motions occur in the transverse nonlinear vibrations of the string-beam coupled system under certain conditions. An erratum to this article is available at .     

Microstructure and Magnetic Properties of Carbon‐Coated Nanoparticles

Abstract

In the present work, microstructure and superparamagnetic properties of two types of carbon‐coated magnetic Ni and Fe nanoparticles [Ni(C) and Fe(C)] are reviewed. High‐resolution transmission electron microscopy (HRTEM), electron diffraction (SAED), and x‐ray diffraction (XRD) analyses have been used to reveal the distinct structural morphologies of Ni and Fe nanoparticles. Moreover, novel carbon‐coated Ni nanoparticle assemblies offer us great opportunities for studying the mechanism of superparamagnetism in particle assemblies. Magnetization measurements [M(T) and M(H) curves] for assemblies of Ni nanoparticles indicate that modified superparamagnetic properties at T > T _B, have been found in the assemblies of Ni(C) particles. The blocking temperature, T _B, is determined to be near 115K under a certain applied field. Above T _B, the magnetization M(H, T) can be described by the classical Langevin function L using the relation, M/M _s (T = 0) = coth (μH/kT) ? kT/μH. It is suggested that these assemblies of carbon‐coated Ni nanoparticles have typical single‐domain, field‐dependent superparamagnetic relaxation properties. Finally, Mössbauer spectra and hyperfine magnetic fields at room temperature for the assemblies of Fe(C) nanoparticles confirm their distinct nanophases that were detected by structural analysis. Modified superparamagnetic relaxation is observed in the assemblies of Fe(C) nanoparticles, which is attributed to the nanocrystalline nature of the carbon‐coated nanoparticles.     

基于光晶格的超冷原子量子模拟

文章从超冷原子研究的视角出发,回顾了用"从下到上"的方案来开展量子模拟研究的历史。超冷原子作为宏观量子态,各个自由度精确可控,是量子模拟的绝佳平台。光晶格将冷原子物理和凝聚态物理融合起来,是其中最重要的技术之一,为超冷原子量子模拟提供了一个扎实的落脚点。近年来,关于拓扑量子模拟的研究日益兴起,成为超冷原子量子模拟新的重要方向。文章介绍这方面近期的一些工作进展。最后分享作者对超冷原子量子模拟的一些思考。     

Extended Lattice Boltzmann Model

Conventional lattice Boltzmann models for the simulation of fluid dynamics are restricted by an error in the stress tensor that is negligible only for small flow velocity and at a singular value of the temperature. To that end, we propose a unified formulation that restores Galilean invariance and the isotropy of the stress tensor by introducing an extended equilibrium. This modification extends lattice Boltzmann models to simulations with higher values of the flow velocity and can be used at temperatures that are higher than the lattice reference temperature, which enhances computational efficiency by decreasing the number of required time steps. Furthermore, the extended model also remains valid for stretched lattices, which are useful when flow gradients are predominant in one direction. The model is validated by simulations of two- and three-dimensional benchmark problems, including the double shear layer flow, the decay of homogeneous isotropic turbulence, the laminar boundary layer over a flat plate and the turbulent channel flow.